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Abstract 

Nucleon structure functions are shown to have a qualitative (or 'formal') relation 
to the classical elliptic theta functions. In particular, 9[/2 shows a clear resem- 
blance to a non-singlet structure function like xF% . In the appropriate range, the 
Q 2 -dependence of the moments of 9[/2 is in near-quantitative agreement with 
QCD, and at low-Q 2 the moments converge to a common value, as observed 
empirically for xF%. At very high Q 2 (Q — ► Mpl), #i/2 — > rf(x) (where r](x) is 
the Dedekind eta function) while xF^ in the same limit appears closer to rj 5 (x). 
A comparison of the theta function identity 0f (0) = 0f (0) - 0f (0) with the rela- 
tion xF% = Q(x) — Q(x) suggests that singlet structure functions have more in 
common with 63 and 84. The possibility of interesting large-x/small-x 'duality' 
relations for structure functions emerges naturally from the analysis. 



1 E-mail:w. g.scott@rl.ac.uk 



1 



Structure functions have long been recognised as fundamental by experimentalists, 
phenomenologists and theorists alike, and impressive progress has been made over 
the years in achieving a precise, quantitative understanding of their behaviour, over 
an increasingly large a kinematic range, mainly within the context of perturbative 
QCD (see for example Ref. |[J). A quantitative description/representation of structure 
functions to match that achieved to date in perturbative QCD, however, is not the 
immediate (nor indeed the ultimate) aim of the present paper, and neither in fact is 
QCD our point of departure here. Rather the aim of this paper is to draw attention 
to a rather startling (if qualitative) empirically-observed similarity /relation existing 
between nucleon structure functions and the classical elliptic (or Jacobi 0) theta 
functions, almost certainly involving, it would appear, the very highest energy scales, 
eg. the unification scale or the Planck-energy. 

Theta functions of various kinds, and in particular the modular symmetry inher- 
ent to theta functions, are themselves recognised as fundamental, mainly within the 
context of more-abstract theory, most notably in string-theory ||. Triggering the 
present paper, and drawing strongly on considerations of modular symmetry (with 
specific reference to the Jacobi theta functions in the scale-invariant case) are recent 
results M on strong/weak-coupling duality, monopole condensation etc, which are gen- 
erally believed to bear on the confinement problem, and which in the long-term might 
reasonably be expected to lead to practical calculations of relevant non-perturbative 
quantities in theories akin to QCD. Thus while (as we shall see) there do seem to be 
points of contact from our observations to perturbative QCD and to a number of asso- 
ciated phenomeno logical ideas/models, we do have to draw attention to the possibility 
of an important link to what one might call 'fundamental theory' here. Naturally, we 
continue to emphasise the empirical aspects of our observations, which we regard as 
striking and interesting in their own right, regardless of explanation. 

With duality ideas in mind therefore, we begin by observing that the modular 
transformation formula (Weyl inversion) for a simple theta function 9[t] viz: 

e\r] = f-e[-i/r] (1) 

is arguably just the kind of relation one might hope to establish between the strong- 
coupling (r — > 0) and weak-coupling (r — > oo) domains in a theory like QCD. The 
parameter r is in general complex (with Im r > 0) but for definiteness we focus here 
on the case that r is purely imaginary: 

r = it (t > 0) (2) 

so that t = —ir is real and positive in this case (via the notion of the 'running' 
coupling, we might think of t as the logarithm of some 'energy', see below). 
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The four classical theta functions 6i(z\r), i = 1—4 |2j all have well-defined modular 
transformation properties, each satisfying a relation similar to Eq. 1 in the case that 
z = 0. In Ref. the functions 62, 63 and 84 (evaluated at z = 0) arise in connection 
with the analysis of spin structures in massless S-dual theories with dimensionless 
couplings. In this paper we make particular use of the function 9[(z\r): 

00 

B' x {z\t) =2^(-l)V 27rT(2n+1)2/8 (2n + l)cos(2n + l)2 (3) 

n=0 

obtained from 81 by differentiation with respect to z (the Qi are by no means all 
independent, eg. 9 2 (z) = e^z + n/2), 6 3 {z) = e^z + ix/2) and 6[(0) = # 2 (0) # 3 (0) # 4 (0) 
H, where &i(0) denotes a theta function evaluated at z = 0). 

With regard to the experimental data on structure functions, theoretical fits based 
on perturbative QCD (which makes rather clear-cut predictions for the Q 2 dependence 
of structure functions, at sufficiently high Q 2 ) have typically been used to interpolate 
and extrapolate existing experimental data to cover the full kinematic range of cur- 
rent/future experimental interest (ie. up to and including LHC energies). By way of 
example, in Figure la we show (solid curves) the x-dependence of the isospin-averaged, 
non-singlet neutrino-nucleon structure function xF%, based on the MRS || fit, for 
factor-of-ten increments in Q 2 , for Q 2 = 1 — 10 7 GeV 2 (here Q 2 is the square of the 
four- momentum transfer and x is the usual Bjorken scaling variable). The data used 
in these fits cover the range Q 2 = 1 — 10 4 GeV 2 , the fit being deliberately restricted to 
Q 2 > 1 GeV 2 , to reduce the risk of trespass into the non-perturbative regime. While 
the x-dependence at fixed Q 2 is largely determined by the data in these fits (pertur- 
bative QCD has much less to say about the expected x-dependence) , a number of 
phenomenological ideas/models do serve as some guide in selecting appropriate func- 
tional forms for fitting. Regge-theory ||, for example, suggests an x^ dependence for 
xF 3 at small x and dimensional counting rules 0] suggest (1 — x) 3 for the behaviour 
as x — > 1, these predictions being expected to hold at some not-too-large Q 2 (Q 2 ~ 1 
GeV 2 ). The sequence of curves shown in Figure la prompts the question of what is 
the high energy limit of the structure function or (perhaps more relevantly) what is 
the extrapolated form at the unification scale or the Planck energy? 

The undenied continued sucess of thermodynamic/statistical-models H f| of struc- 
ture functions would certainly suggest that statistical-mechanical (or combinatorial) 
considerations could underlie the form of the structure functions at high energies (it 
has also been noted directly from experiment that the x-dependence of structure func- 
tions becomes more exponential in character with increasing energy JCTI). Closely 
associated with theta functions, the Dedekind eta function 7] [t] obeys the same mod- 
ular transformation formula as the simple theta function (Eq. 1 above) and has a 
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clear combinatorial significance (see below). Furthermore the x- dependence of the eta 
function, given by: 

i 00 

■q{x) = £24 Y[ (1 - x n ) (4) 

n=l 

(t](x) = r/[r], a; = exp(27ur)) is remarkably reminiscent of the x-dependence expected 
of a non-singlet structure function, cf. the Regge prediction etc. (above). With t = —it 
taken to be real and positive as before (cf. Eq. 2) we have: 

t = —In- (5) 

with x also real and positive and satisfying < x < 1 as required. With this (perhaps 
unexpected) association of t with the scaled parton energy x (rather than Q = y/Q 2 ), 
the modular transformation (cf. Eq. 1) would become a large-x/small-a; relation, ie. a 
symmetry versus hit = ln^-lnl/x centered about a 'critical' x- value xq = exp(— 2n). 

In Eq. 4 the pre-factor x 1 / 24 is known to be essential to the modular covariance. 
The infinite product has a clear combinatorial significance as the reciprocal of the 
generating function for the classical partition function [IT]. Physically it is just the 



reciprocal of the ordinary partition function for a system comprising a (freely) variable 
number of indistinguishable massless bosons moving freely in one dimension (the index 
1/24 in Eq. 4 is then interpreted as the associated Casimir (or vacuum) energy WW). 
The generalisation to more than one 'colour' of boson is readily obtained by raising the 



generating function to the power of the number of boson colours ||13|| . Functions like 
7] 3 (x), rj s (x) . . .appear naturally as operator traces or 'characters' in relation to the 
infinite dimensional lie algebras su 2 , su 3 ...which are equivalent to (2-dimensional) 
current algebras The (one-colour) generating function for fermionic partitions, 

related to r](x)/r](x 2 ), is the same infinite product (Eq. 4) but with {l+x n )~ l replacing 
(1 — x n ). The bosonic entropy is then \/2 times the fermionic entropy and (all else 



being equal) bosons take twice the energy of fermions aymptotically |TJ 

In a heuristic spirit, we are in effect arguing that the above functions relate to 
the probability (per unit logx) of finding the system in a particular state, and as 
such represent structure functions in the high energy limit, if we assume that partons 
behave like free particles at high energy. In Figure 2a we plot the moments of xF 3 (solid 
circles) for factor-of-ten increments in Q 2 (as above) calculated from the QCD fit, as a 
function of lnQ 2 /A 2 on a log-log plot, ie. we plot lnM 3 (iV) versus s = lnlnQ 2 /A 2 (with 
A ~ 250 MeV) so that perturbative QCD evolution is reduced to a linear dependence 



with a known slope (M 3 (N) = J x N l F%dx [I5|). Extrapolating (the moments of) 



xF 3 to Q 2 ~ 10 38 GeV 2 (Q ~ M PL ), we found to reasonable accuracy: xF 3 — > 
[xi2r](x)/r)(x 2 )} 5 , where the moments of this function are represented in Figure 2a by 
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the open half-circles plotted at the largest value of Q 2 . This function (clearly cloaely 
related to the fermionic partition function above) is plotted versus x in Figure la 
(dashed curve, labelled Q 2 ~ M| L ) and is in fact very similar in general appearance 
to the function rf(x) (not shown), either function satisfying the GLS sum-rule [IB] for 
the N = 1 moment to an accuracy of a few percent. 

The above observations (Figure la/2a) lend empirical support to the notion that 
powers (or combinations of powers) of 77 (x) may underlie the functional form of the 
structure functions at the highest energies (regarding the basic shape and also, it 
would seem in this case, the normalisation). On the other hand, we have no particular 
explanation for the apparent preference for the fifth power here (but see Ref. |31 
below) . 

The theta function 6[ (Eq. 3) evaluated at z = is proportional to the third power 
of the eta function, specifically: 

9' 1 (0,x)/2 = r ] 3 (x). (6) 

where we adopt the usual convention 9i(z,x) = 9i(z\r) (the conventional symbol q 2 
often used in place of x, is deliberately avoided here). In Figure lb we plot the 
evolution of 6[ for factor-of-ten increments in Q 2 , with the a posteriori (see below) 
identification (assumed valid in the limit z — > 0): 

z ~ const. - ^lnlnQ 2 /A 2 ~ 0.6) (7) 

Z7T 

where the constant is fixed such that z = corresponds to Q 2 ~ 10 38 GeV 2 (ie. the 
constant is equal to the Q 2 -dependent term with Q set equal to Mpl, viz. const. ~ 
(l/27r) In In Mp L / A 2 where M PL is the Planck-mass, M PL ~ 10 19 GeV). Comparing 
Figure la and Figure lb, it is apparent that the ^-dependence of the theta func- 
tion shows a really rather startling similarity to the Q 2 dependence of the structure 
function, at least over the current experimentally relevant range Q 2 = 1 — 10 7 GeV 2 
{z - 0.55 - 0.25). 

A more quantitative comparison is achieved by plotting (the logarithm of) the 
x-moments of the theta function versus z (Figure 2b) alongside the corresponding 
plot (Figure 2a) for the structure function moments versus s = lnlnQ 2 /A 2 . In each 
case we see a broadly linear dependence of (the logarithm of) the moments on the 
ordinate, with a slope which increases with increasing moment number. The theta 
function moments (Figure 2b) converge to a common value (2ir) as z — > n/2, where 
6[/2 — > 2tt5(x — 1). While the structure function moments cannot be plotted (versus 
s) for Q 2 < A 2 and are not plotted anyway for Q 2 < 1 GeV 2 , as explained above, it 
is evident nonetheless (Figure 2a) that if the linear dependence of the QCD fit were 
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naively continued to lower Q 2 (broken lines) the structure function moments would 
also apparently converge in very similar way (this striking feature of the data has no 
explanation in the context of perturbative QCD, and has been cited as evidence for 
'valons' ]T7| , and is also apparently a prediction of bag-models |TB| and possibly also 
chiral-soliton models |0|)- Of course the structure function becomes proportional to 
5(x — 1) at very small values of Q 2 , as a consequence of the dominance of elastic 
scattering (xF 3 -> 5.7 x 8(x - 1) , as Q 2 -> []T3J). 

Assuming the correspondence Eq. 7, in Figure 3 we compare the slopes (with 
respect to logQ 2 ) of the moments of the theta function (solid curve) directly with 
the leading-order QCD non-singlet anomalous dimensions, evaluated for three quark 
flavours (broken histogram) for the first ten moments (N = 1 — 10). The theta function 
slopes are evaluated at some not-too-small value of z (z = 7r/6 ~ 0.5), roughly where 
present-day measurements cluster. Somewhat remarkably the theta function slopes 
are seen to reproduce the three flavour anomalous dimensions to an accuracy of a few 
percent or so, over the range of moment numbers plotted, and it is this comparison 
which justifies, (a posteriori) the identification Eq. 7, where the normalising factor 
l/(2ir) is inserted purely empirically. At large N the non-singlet anomalous dimensions 
grow like ||logAf, while the corresponding dependence for the theta function slopes 
appears to be -y/2N. Large N moments are not explored experimentally, however. 

Note that this remarkable quantitative agreement, whether partly (or even wholly) 
'accidental', extends over the entire kinematic range explored/fitted experimentally. 
At higher values of Q 2 (z < 7r/12) the z-dependence of the theta function moments 
(Figure 2b) begins to 'flatten-out', with Bjorken scaling holding 'exactly' for the theta 
function at the very highest energies (Q —>■ Mpi). While this effect is naturally not 
present in the QCD fit (extrapolated to high energy) it is clear that such high energies 
(Q 2 > 10 7 GeV 2 ) have so far also not been explored experimentally. (Certainly in any 
case, QCD can hardly be expected to hold unmodified for Q — ► M PL ). Taking the 
Q 2 -dependence of the theta function seriously, at the lower values of Q 2 (z > tt/12) 
these effects are small (see above) and better than 1% measurements of xF 3 would be 
required to establish them even at HERA, Q 2 ~ 10 4 GeV 2 [pfj|| . By the same token, 
substantive deviations (^ 20%) from the familiar perturbative QCD Q 2 -dependence 
might be expected for Q 2 > 10 7 GeV 2 and could conceivably show-up, for example at 
a next-generation ep collider. 

To summarise, we have seen that the classical theta function 8[ bears a clear 
resemblance to a non-singlet structure function, like xF 3 . While these two functions are 
not identical, even in the region explored experimentally, it is hard to believe that the 
observed similarities are entirely accidental. We are led to conclude that a qualitative 
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relationship exists between structure functions and theta functions, a notion on which 
we will enlarge below. This is the main conclusion of our investigation and it rests on 
the analysis of the non-singlet structure functions, Figures 1-3 above. 

Regarding singlet structure functions, ie. F 2 = Q(x) + Q(x), or the closely related 
gluon distribution G(x) in QCD, we have only a few comments. First, we note a 
similarity of form between, for example, the relation xF% = Q(x) — Q(x) and a theta 
function identity like 0%(O) = 0f(O) — 0f (0), essentially the Jacobi 'obscure' formula 0: 

oo oo oo 

I6x? n (! + xn T =11(1 + x n ^f - A (1 - x n -^ f (8) 

n=l n=l n=l 

Multiplying the theta function identity by the infinite product ]\{l — x n ) 2A and plotting 
the x-dependence of the various terms in Figure 4, the resulting curves are reminiscent 
of the familiar behaviour of Q(x), Q(x) etc. at small x. This suggests that singlet 
structure functions may be related in some way to functions like 83 and 84 which are 
expressible as simple infinite products @, without any powers of x. In this connection, 
we note further that, exploiting the modular covarince of Eq. 4, the infinite product 



pl(l — x n ) can be shown [] 1] to be remarkably well approximated by: 

^2 



v^-^= T exp(-^- T ) (9) 



X 



(with less than 10% deviation even for X cLS small cLS X — 10 6 for example) where the 
right-hand side bears a clear resemblance to the well-known Q 2 — > 00 expression for 
the (unintegrated) gluon distribution at small x (x < 1CT 3 ) pl| : 



dG(x) x ->o 



7 exp(- -g 7 ^fc T ) (10) 



d\nk 2 T \kl J^J PV 56C(3)alni 



obtained by solving the BFKL |22|] equation analytically |23j (and it will be noted that 



despite the 'formal' inconsistency of the limits (Eq. 9 vs. Eq. 10) the ranges of validity 
actually overlap over three orders of magnitude in x). In Eq. 10, a = (where a 
is the QCD gauge- coupling) and Xbfkl = 4 In 2 a (k 2 ,, k\ are non-perturbative scale 
parameters). 

It should be pointed out that if the analytic similarity of Eq. 9 to Eq. 10 is not 
accidental, then there may be someting to learn here about how the BFKL singularity 
is regulated in QCD (Eq. 9 LHS does not obviously violate the Froissart bound). 
Likewise, the gauge-coupling |23] would seem to be related to the index 1/24 above (we 



note that for QCD, the Freudenthal de Vries 'strange' formula reduces to 8 x Ca = 24, 
where Ca = 3, while for three generations of quarks the corresponding net strength of 
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the quark-gluon vertex in QCD is given by 18 x Cp = 24, where Cp = 4/3). Casimir 
himself postulated a relation between the zero-point energy coefficient and the gauge 
coupling P3| . In (2-dimesional) dilaton gravity, the conformal anomaly c/24 is known 
to play the role of the (dimensionless) inverse Newton constant ||26|| . 

Of course theta functions solve the diffusion equation d 2 9/dz 2 = (Ai/tt) dO/dr and 
the link between BFKL and diffusion has long been recognised, having its physical 
origin in the absence of transverse momentum (kp) ordering in the BFKL ladder [[27[ . 
While the kind of global similarity to theta functions seen in the non-singlet case 
(Fig. 1-2) would certainly suggest that diffusion plays an important role for longitu- 
dial variables also, clearly it is not possible (from perturbative QCD alone) to derive a 
double-differential equation for structure functions (comparable to the diffusion equa- 
tion) for arbitrary x and Q 2 . In the limit Q 2 — > oo, x — > the DGLAP equation 
for G(x) in fact reduces to a wave equation: d 2 G/dzdr ~ 2ndG/dr + . . but 
it should be remembered that the DGLAP approach omits explicitly the unordered 
hp contribution above (assumes strong kp ordering). Regarding our extrapolation to 
Q — > Mpi, where presumably quantum gravity effects cannot be neglected (ie. beyond 
the QCD context) it is perhaps worth noting that the partition function for a micro- 
scopic black-hole is known to satisfy the ID-diffusion equation While the smooth 
extrapolation of DIS physics to black-hole physics implied in the present paper should 
obviously be viewed with due caution, we note further that a partial analogy between 
microscopic black- hole evaporation (via Hawking radiation [|3T|) and QCD fragmenta- 



tion has in fact already been recognised in the literature |32| (and see also Ref. |34 
below) . 

If anything, the above observations on singlet structure functions reinforce the 
case for a connection between structure functions and theta functions, which was 
the main conclusion of the non-singlet analysis. That is not to say that we have 
in hand a complete and quantitative description of the neutrino-nucleon structure 
functions in terms of the Jacobi theta functions. The comparisons we have made are 
at best only partially successful and sometimes apparently self-contradictory (eg. xF 3 
cannot 'equal' 9[ which behaves like x^ at small x and simultaneously display an x^ 
dependence as implied by Eq 8). Also some generalisation of Eq. 7, to cover the low Q 2 
region, needs to be found, which may be difficult in practice. Nonetheless, considering 
the non-singlet and singlet cases together, the accumulated evidence for a qualitative 
(or 'formal') relation between structure functions and theta functions is undeniably 
rather strong. The most obvious possibility would be that structure functions are theta 
functions of some sort (or powers of theta functions) or some variant of theta functions 
with similar behaviour. Possibly, theta functions relate to some 'idealised' structure 
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function eg. of a pure-gauge system (or in the massless/SUSY limit) unobservable 
experimentally. The theoretical derivation of such a result might conceivably proceed 
within the representation theory of some infinite dimensional algebra ]T3| . 

To conclude, our point of view is that theta functions arise in connection with the 
highest energy scales, eg. the Planck-scale or the unification scale (if our observations 
mean anything at all, some high energy scale must be associated with the z = for 
the theta functions, and the Planck/unification scale seems the obvious best guess, 
cf. Eq. 7). While in the QCD context we are often used to evolving models/fits from 
lower to higher energies, one might reasonably argue that symmetry and simplicity are 
in fact more naturally associated with the highest energy scales, so that ultimately, 
evolving functional forms down in energy may come to seem much more sensible. We 
remark once more that modular functions and modular symmetry already feature 
prominently in theoretical considerations of Planck-scale physics. For example, string 
theory threshold corrections to the reciprocal gauge-couplings |33| turn out to be in 



part proportional to ln|^(r)| 4 , (with r the value of a given moduli field scaled to 
the Planck-energy). Also, in five- dimensional D-brane theory the partition function 
of a microscopic black-hole (small with respect to the compact ification radii) can be 
expressed in terms of integral powers of the eta-function [33]. The possibility that 



such apparently abstract theoretical concerns may already be mirrored in existing 
phenomenology (cf. xF 3 — > ~?7 5 (a;) for Q — > M PL ) is both exciting and potentially 
important. 

For the future, emerging naturally from this analysis large- a;/small-x 'duality' re- 
lations for structure functions of the form Eq. 1 (or similar) are certainly an attractive 
and ultimately testable possibility. While the non-singlet case would seem to offer the 
best chance, empirical fits to singlet structure functions (eg. to F 2 ep Ref. do in fact 
suggest 'critical' x-values of a sort, apparently Q 2 -independent (with x ~ exp(— n)). 
With some correlation between the large- a; and small- x behaviour of structure func- 
tions undoubtedly physically plausible |36[ , we have every reason to hope that definitive 



experimental or theoretical evidence for the existence of relations of the above form 
will eventually be found. In this last regard, we are pleased to call attention to the 
work of R. Janik |I7|, pointing to a connection between xF 3 and the classical elliptic 
theta functions, established (for transverse variables) in the Regge limit of QCD, via 
the modular invariance of the 'odderon'. 
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Figure Captions 



Figure 1 a) The x-dependence of the non-singlet neutrino-nucleon structure function 
xF 3 , based on the standard MRS QCD fit to existing data, for factor-of-ten incre- 
ments in Q 2 , Q 2 = 1 — 10 7 GeV 2 (solid curves). The dashed curve labelled Q 2 ~ M| L 
is the function [x^r/^x)/^^ 2 )} 5 (where T](x) is the Dedekind eta function [pl| Eq. 4) 
which provides a reasonably accurate representation of the QCD fit extrapolated to 
Q 2 ~ 10 38 GeV 2 (see Figure 2a). b) The theta function 9[(z, x)/2 plotted as a function 
of x (solid curves) for a sequence of z- values, z = 0.55 — 0.25 (see Eq. 7). The limiting 
form of 6[/2 at z = is rf(x) (dashed curve). 

Figure 2 a) (The natural logarithm of) the x-moments (N = 1 — 7) of the structure 
function xF 3 , plotted as a function of s = lnlnQ 2 /A 2 , for Q 2 = 1 - 10 7 GeV 2 (filled 
points). Extrapolating the QCD dependence (solid lines) to Q 2 ~ 10 38 GeV 2 , the 
extrapolated moments are seen to be reasonably accurately reproduced by the mo- 
ments of the function [x^r](x)/r](x 2 )] 5 (open half-points). A naive extrapolation to 
low Q 2 (broken lines) shows the moments apparently converging (see text), b) The 
^-moments of the theta function (9' 1 /2) plotted as a function of z. The theta function 
moments also converge to a common value as z — ► 7r/2 (ie. at small values of Q 2 ) 
where 0[/2 — > 2tt5(x — 1). At small z (ie. at large values of Q 2 ) the z-dependence 
'flattens-off, so that Bjorken scaling holds exactly for 8[ as z — > 0. 

Figure 3 The slopes of the theta function moments with respect to logQ 2 , (assuming 
the correspondence Eq. 7) evaluated at z = n/6 ~ 0.5, plotted as a function of 
moment number (solid curve). The broken histogram shows the leading-order QCD 
anomalous dimensions, calculated for three quark flavours. The theta function slopes 
reproduce the leading-order QCD anomalous dimensions (to an accuracy of a few 
percent or so) over the range of moment numbers plotted. 

Figure 4 Terms in the theta function identity #|(0) = ^1(0) ~ 9f(0) @ (essentially the 
Jacobi 'obscure' formula Eq. 8) multiplied by the infinite product 11(1 — x™) 24 , plotted 
as a function of x, for small-x. The resulting curves are reminiscent of the familiar 
behaviour of Q(x), Q(x) etc. at small-x (cf. # 4 (0) = 0f (O)-0|(O) : xF 3 = Q(x)-Q(x)), 
suggesting that singlet structure functions may be related to functions like 6*3 and #4, 
expressible as simple infinite products, without any powers of x. 



13 




II 



CO 

II 




Figure 5 




0.05 0.10 0.15 



Figure 4 



